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Abstract 

We obtain the affine Euler-Poincare equations by standard Lagrangian reduction 
and deduce the associated Clebsch-constrained variational principle. These results 
are illustrated in deriving the equations of motion for continuum spin systems and 
Kirchhoff 's rod, where they provide a unified geometric interpretation. 
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1 Introduction 

Lagrangian reduction by symmetry (LRS) applies when Hamilton's principle is invariant un- 
der the action of a Lie group. The implications of Lie symmetries of Hamilton's principle sum- 
mons many mathematical concepts in a wealth of interesting applications of geometric me- 
chanics in the dynamics of particles and fields, as well as in control theory For example, the 
Euler-Poincare (EP) reduction framework discussed in Cendra, Holm, Marsden, and Ratiu 
[1998] and Holm, Marsden and Ratiu [1998] summons all of the adjoint and coadjoint ac- 
tions of Lie groups on themselves, both on their representation vector spaces and on their 
Lie algebras. Several applications of the EP reduction framework to classical continua are 
surveyed in Holm, Marsden and Ratiu [1998]. 

Continuum condensed matter theory deals with motions of a wide range of complex 
fluids; that is, fluids with internal degrees of freedom, e.g., phases, spins, directors or other 
order parameters. Order parameters are geometrically either objects in a vector space, or 
coset spaces of broken Lie symmetry groups with respect to subgroups that leave these 
objects invariant. The evolution of the order parameters of a complex fluid is called its 
micro-motion. Complex fluids whose order parameters are continuous material variables are 
said to be perfect complex fluids (PCFs). In contrast, discontinuities and singularities in 
the order parameter of a complex fluid are called defects and these may either be frozen 
into the fluid motion, or they may move relative to the material. The classical examples 
of complex fluids are liquid crystals, magnetic materials, and the various types of Landau 
theories of superfluids. The application of the EP reduction framework to the classical PCF's 
is surveyed in Holm [2002]. 

An interesting feature arises for the EP dynamical formulation of a perfect complex fluid 
whose order parameter is carried along by the flow of the material parcels. Namely, the 
Legendre transformation to its Hamiltonian formulation yields a noncanonical Lie-Poisson 
bracket that is dual to a Lie algebra possessing generalized 2-cocycles. The generalized 2- 
cocycles appear in the corresponding Lie-Poisson Hamiltonian operator as covariant deriva- 
tives with respect to a Lie-algebra valued connection, in precisely the same form as found 
in the Yang-Mills fluid plasma for chromohydrodynamics Gibbons, Holm and Kupershmidt 
[1982, 1983], in superfluid Helium with vortices Holm and Kupershmidt [1982, 1987], in 
type-II superconductors Holm and Kupershmidt [1983a] and in the classical dynamics of a 
spin glass Holm and Kupershmidt [1988]. Following Holm and Kupershmidt [1983b] and 
Marsden and Weinstein [1983], the same Lie-Poisson Hamiltonian operator was also obtained 
in Holm [2002] by applying Clebsch constraints to Hamilton's principle to enforce the aux- 
iliary kinematic equations obeyed by the PCFdynamics. While these previous approaches 
revealed the presence of the generalized 2-cocycles, they did not specifically reveal their 
source. However, recent developments in Gay-Balmaz and Ratiu [2009] provide the means 
of deriving these generalized 2-cocycles by showing that they are summoned when the Lie 
symmetry of Hamilton's principle is an affine action. Thus, affine Lie symmetry reveals the 
source and guarantees the presence of generalized 2-cocycles. The present paper pursues 
these matters further for PCF's and for complex fluids whose defects are frozen into the 
flow. 

Reduction of Hamilton's principle by affine Lie symmetry is used here to construct 
a mathematical framework for deriving the affine Euler-Poincare equations governing the 
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nondissipative coupled motion and micromotion of complex fluids. More precisely, we show 
how the affine Euler-Poincare reduction theorem can be seen as a particular case of the gen- 
eral Lagrangian reduction process TQ — > TQ/G, that is, the case of a Lagrangian defined on 
a tangent bundle, and invariant under the tangent lift of a Lie group action. This is achieved 
by modifying the Lagrangian, in the same spirit as it is done for the Euler-Poincare reduc- 
tion in Cendra, Holm, Marsden, and Ratiu [1998]. As a consequence, we obtain and explain 
the associated Clebsch-constrained variational principle in a natural way. As examples of 
physical interest, we apply the new variational principles to continuum spin systems and to 
the Kirchhoff rod. 

The plan of the paper is as follows. In the first part of Section 2 we recall the relevant 
facts about affine Euler-Poincare reduction as it applies to complex fluids and review the 
process of general Lagrangian reduction. Then it is shown how to obtain the affine Euler- 
Poincare equations by standard Lagrangian reduction and the associated Clebsch-constrained 
variational principle is presented. In Section 3, we apply the resulting theory to the dynamics 
of continuum spin systems. In particular, we get the variational principle 

5 J (Z(i/, 7 )+w(7-cPV)) = 0, 

for variations 8j, 5w, and constrained variations 5v = rj+ [u, 77]. In Section 4, we review some 
needed facts from Kirchhoff 's theory of rods, using the Kirchhoff and spatial representations, 
and we show how the body representation can be obtained by affine Euler-Poincare reduction, 
in the case when central potential forces are considered. Finally, by applying the theory 
developed in Section 2 we obtain the Clebsch-constrained variational principle for Kirchhoff 's 
equations, 

5 J ^(w,7,n,r,p) + y u- (tt - 8 s uj - n, x w ) 

+w- (r + w x r - d s j- Ct x 7) + f(p + w x p-j) ^ dt = 0, 

for variations 5Q, 5T, Sp, Su, 8w, 5f and constrained variations 

drj dv 
ou) = — — h uj x 77 and 07 = — — \- uj x v — 77x7. 

ot at 



2 The affine Euler-Poincare Equations via Lagrangian 
reduction 

The goal of this section is to show how the affine Euler-Poincare principle can be obtained 
by standard Lagrange reduction of a cotangent bundle. From this alternative approach to 
the affine Euler-Poincare equations we obtain, in a natural way, the associated Clebsch- 
constrained variational principle. 

We begin by recalling some needed facts about affine Euler-Poincare reduction for semidi- 
rect products (see Gay-Balmaz and Ratiu [2009]). Let V be a vector space and assume that 
the Lie group G acts on the left by linear maps (and hence G also acts on the left on its dual 
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space V*). As a set, the semi direct product S = G(§)V is the Cartesian product S = G x V 
whose group multiplication is given by 

(gi, U1X02, v 2 ) = (gig 2 , v x + gxv 2 ), 

where the action of g G G on v G V is denoted simply as gv. The Lie algebra of S is the 
semidirect product Lie algebra, S = Q®V, whose bracket has the expression 

ad( &>Wl )(6,V2) = (£2,^2)] = ([&, 6]> -6^1), 

where £v denotes the induced action of g on V, that is, 



iV '-=dt 



exp(t£)v G V. 



t=o 



From the expression for the Lie bracket, it follows that for (£, i>) G s and (/x, a) G s* we have 

adj^^/x, a) = (adg fi — v o a, — £a), 
where y* and 1100G5* are given by 



d 
dl 



exp(t£)a and (voa,£) l 



■(£a,v) Vl 



t=o 



and where (-, •) : g* x g — > K. and (•, -) y : V* x V — > K are the duality parings. 

Given a Ze/t representation of G on the vector space V*, we can form an affine left 
representation 9 g (a) := ga + c(g), where c G .F(G, V*) is a left group one-cocycle, that is, it 
verifies the property 

c(gh) = c(g) + gc(h) (2.1) 

for all g, h G G. Note that 

#cx P (to( a ) = £ a + dc (0 



t=o 



and 



(£a + dc(f), v) v = (dc T (v) - v o a, £) , 
where dc : g — * V* is defined by dc(£) := T e c(£), and dc T : — > g* is defined by 

(dc T (v),O s --=(dc(0,v) v . (2.2) 

Affine Lagrangian semidirect product theory. 

• Assume that we have a function L : TG x V* — > K. which is /e/t G-invariant under the 
affine action (v h ,a) i-> (gv h ,6 g (a)) = (gv h , ga + c(g)). 



In particular, if ao G V"*, define the Lagrangian L ao : TG — > K. by 

L ao (v g ) := L(Vg,a ). 



(2.3) 



Then L ao is left invariant under the lift to TG of the left action of G c ao on G, where 
Gjj o is the isotropy group of qq with respect to the affine action 9. 
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• Left G-invariance of L permits us to define / : g x V* — > M by 

l(g~~ x v g ,6 g -i{a Q )) = L(v g ,a ). (2.4) 

• For a curve g(t) G G, let £(£) := g{t)~ l g(t) and define the curve a(t) as the unique 
solution of the following affine differential equation with time dependent coefficients 

d = -£a-dc(0, (2.5) 
with initial condition a(0) = ao- The solution can be written as a(t) = 6 g ^-i(ao). 

Theorem 2.1 With the preceding notations, the following are equivalent: 

i With ao held fixed, Hamilton's variational principle 

5 ^ L ao (g,g)dt = 0, (2.6) 

J to 

holds, for variations 8g(t) of g(t) vanishing at the endpoints. 

ii g(t) satisfies the Euler- Lagrange equations for L ao on G. 

iii The constrained variational principle 

5 [ 1 Z(f, a)dt = 0, (2.7) 

/io/ds on q x V* , upon using variations of the form 

dri 

6£ = — + [£, r]\, 5a = -rja - dc(r/), 

where r/(t) G Q vanishes at the endpoints. 

iv The affine Euler-Poincare equations hold on g x V* : 

d 61 „6l 61 rp / 6l\ , no . 

adj — + — o a - dc 1 — . (2.8) 



dt6i * 8£ 5a \5a 

See Gay-Balmaz and Ratiu [2009] for a proof and applications to complex fluids. The 
affine Euler-Poincare equation (2.8) can be rewritten in the form of a conservation law as 
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General Lagrangian reduction. Consider a left action $ : GxQ- > Q of a Lie group G 
on a manifold Q. Let L : TQ — > R be a G-invariant Lagrangian under the cotangent-lifted 
action of G on TQ. Because of this invariance, we get a well defined reduced Lagrangian 
I : (TQ)/G -> R satisfying 

/(K])=L(^). 

Assuming the group action is free and proper, the quotient space (TQ)/G is intrinsically a 
vector bundle over T(Q/G) with a fiber modeled on the Lie algebra g. Using a connection 
A on the principal bundle n : Q — > Q/G we have a vector bundle isomorphism 



a A : (TQ)/G — T(Q/G) © g, [u,] .— au(KD := (Ttt«), 



G 



over Q/G, where the associated bundle g := Q g, is defined as the quotient space of 
Q x g relative to the left action (q,£) Ad 9 £) of G. The elements of g are denoted 

by v = [q, £\q. Using the isomorphism a A , we can consider I as a function defined on 
T(Q/G) © g, and we write l(x,x,v) to emphasize the dependence of I on (x,x) G T(Q/G) 
and v6g. However one should keep in mind that x, x, and v cannot be considered as being 
independent variables unless T(Q/G) and g are trivial bundles. 
We now formulate the Lagrangian reduction theorem. 

Theorem 2.2 The following conditions are equivalent: 
i Hamilton's variational principle 



5 [ L(q,q)dt = 0, 
Jtn 



holds, for variations 5q(t) vanishing at the endpoints. 

ii The curve q(t) satisfies the Euler- Lagrange equations for L on TQ. 

iii The reduced variational principle 

5 / l(x, x, v)dt = 
J t 

holds, for variations Sx © 5 A v of the curve x(t) ®v(t), where 5 A v has the form 

S A v = —fj + [v, fj] + B(Sx, x), 

with the boundary conditions 5x(ti) = and fj(U) = 0, for i = 0,1. 

iv The following vertical and horizontal Lagrange-P oincare equations, hold: 

D dl dl . 

— — [x, x, v) = ad s — {x, x, v) 
Dtdv dv 

dl D dl I dl ~ x 

Tx^ *' ^ -Joidi^ *• ^ = W x ' *' ^ hB{x) 
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We now comment on the various expressions appearing in parts ii and iii. In the expres- 
sion 

5 v = —fj + [v, fj] + B(5x, x), 

D/Dt denotes the covariant time derivative of the curve f](t) e g associated to the principal 
connection A, that is, for fj{t) = [q(t), £(t)]c, we have 



§- t [q(t),m] G = [q(t), i(t) - [A(q(t),q(t)),m 



G 



The bracket [v, ff\ denotes the Lie bracket induced by g on each fiber of g. The two-form 
B G Q 2 (Q/G, g) is the curvature on the base Q/G induced by the curvature form B = 
dA — [A, A] G Q 2 (Q,q) of A. Notice that for the formulation of the Lagrange- Poincare 
equations, the introduction of an arbitrary connection V on the manifold Q/G is needed. 
For simplicity a torsion free connection is chosen. The partial derivatives 

dl dl 

— (x,x,v)e T* (Q/ G) and — (x, x, v) € g* 

ox ov 

are the usual fiber derivatives of / in the vector bundles T(Q/G) and g, and 

dl 

—(x,x,v)eT:(Q/G) 

is the partial covariant derivative of I relative to the given connection V on Q/G and to 
the principal connection A on Q. In (2.10), the covariant derivatives D/Dt are respectively 
associated to the principal connection A and to the affine connection V, on the dual bundles. 
We refer to Cendra, Marsden, and Ratiu [2001] for details and proofs regarding the Lagrange- 
Poincare equations. 



Application to the affine Euler-Poincare equations. Let us assume all the hypotheses 
of Theorem 2.1. We define the Lagrangian 

L : T(G xV*xV)->R, L(v g ,a ,v ,a ,v ) := L(v g ,a ) + (a ,v ), (2.11) 

which is invariant under the tangent lift of the action 

(g, a , v ) i (hg, ha + c(h), hv ) (2.12) 

of G on G x V* x V, given by 

(v g , a , v Q , d , v ) i-> {hv g , ha + c(h), hv , ha , hv ). (2.13) 

The idea in definition (2.11) is to introduce the condition that a® is conserved by making it 
the momentum conjugate to a cyclic variable, just as one does for the charge in Kaluza-Klein 
theory. Notice that the Euler-Lagrange equation associated to the variable vq is 

d = 0. 

Thus <2q is a constant. 
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The variable v$ is not constant, but writing the Euler-Lagrange equation for ao yields 
the first order equation 

r)T 

v -j- = 0. (2.14) 

Thus, the Euler-Lagrange equations for L are equivalent to the Euler-Lagrange equations 
for L ao with the parameter ao fixed (together with the uo-equation (2.14)). 

In order to apply the Lagrangian reduction process (Theorem 2.2) to the Lagrangian L, 
we need to consider the principal bundle 

G x V* x V -> (G x V* x V)/G 

relative to the G-action (2.12). This bundle turns out to be isomorphic to the trivial bundle 

G x V* x V -»• V* x V. 

To see this note that the diffeomorphism 

if>:GxV* xV ->GxV* xV, ^(#,a ,v ) := (g, g^ao + c(g~ l ), g~ l v Q ) =■ (g,a,v), 

is G-equivariant relative to the action (2.12) and left translation of G on the first factor only, 
that is, 

ip(hg, ha>o + c(h), hv ) = (hg, a, v). 

Using the formula 

Tip~\g, g, a, v, a, v) = (g, g, ga + c(g),gv, g(a + g~ x ga + dc(^ -1 ^)), g(v + g~ l gv)), 

we obtain that the composition L v := L o Tip" 1 is given by 

L V (g, g, a, v, a, v) = L(g, g, ga + c(g)) + <d + g~ l ga + dc{g~ x g), v) . (2.15) 

We shall show that the reduced Euler-Lagrange equations for the Lagrangian L v are the same 
as the affine Euler-Poincare equations for /. To do this, we first calculate the reduced bundle 
T(Q/G) © g, for Q := G x V* x V. Since the G-principal bundle Q := G x V* x V -> V* x V 
is trivial, we have g = V* x V x g, and we can choose the trivial principal connection 

A(g, g, a, v, a, v) = gg' 1 . (2.16) 

Thus we obtain 

TQ/G = T(Q/G) © = x T(V* x V) 3 (£, a, v, a, v), 

[g, g, a, v, a, v] i — > (g^g, a, v, a, v). 
The reduced Lagrangian induced by L v is given by 

l v (£, a, v, d, v) = a) + (a + £a + dc{£), v). (2.17) 

Note that we can write 

l v (f , a, v, a, v) = a) + (6 (a, v), (a, v)) + (J (a, v),£), (2.18) 
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where J : V* x V — > g* is the momentum map of the lift of the action a h- > ga + c(g) to the 
cotangent bundle T*V* ^V* xV, that is, 

J(a,v) = -voa + dc T (v) , (2.19) 

and 9q is the canonical 1-form on T*V*. Recall that ($o(a,v), (b,w)) = (b,v) for v,w G V, 
a,b G V*. We now compute the reduced Euler-Lagrange equations (2.10) associated to l v . 
The horizontal equation for a is given by 



d 5l v 51 



v 



dt 5a 5a 

Using the equalities 



0. 



5l v , 5l v 51 _ 

— = v and — = — £v, 

da da da 

we find that this equation is equivalent to 



C7 

v + £v- — = Q. (2.20) 



5a 



The horizontal equation for v is the equation 



d 5l v 51 



v 



0. 



dt 5v 5v 

which is readily seen to be equivalent to the equation 

d + £a + dc(£) = 0. (2.21) 
Since the curvature of the connection (2.16) vanishes, the vertical equation in (2.10) becomes 

d 5l v (5l v \ 

JtH - ad ^ ■ 



However 



5l v 51 T 
—— = v o a + dc ( v ) 



and thus the vertical equation takes the form 



d 51 

Jts~l 



v o a — v o a + dc T (v) = ad^ (j^j ~~ ac ^(^ ° a ~ dc r (v)) 



or, using equations (2.20) and (2.21) 



d 51 51 . . , , , T ( 51 , 

— -r- - — oo + (!jo a+v o^a + vo dc(£) + dc — - dc 1 (£v) 
dt5t, 5a \5a/ 



ad^ (j^j ~ a( ^l( v ° a ~~ dc T (v)). 



Gay-Balmaz, Holm Sz Ratiu Variational principles, spin systems and Kir chhoff rods 10 



Using the equalities 

(£voa + vo£a + vo dc(f) - dc T (£v), r/) = ((£77 - rj£)a, v) + (£dc(r/) - rjdc(£),v) 

= ([C,v]a,v) + (dc([C,v]),v) 
= (-voa + dc T (v), [£,rj\) 
= -(ad^ (voa-dc T (v)) ,n) , (2.22) 

we finally obtain the affine Euler-Poincare equation 

d 51 „6l 51 tT fM\ 
— — : = ad| — : + — o a - dc — . 
dt6£ ? 5£ 5a \5a J 

The following theorem summaries the results obtained so far. 

Theorem 2.3 Let g(t) be a curve in G, fix an element a e V* , and define the curve 
£(£) := (yf(t)~ 1 5f(t). Consider the left invariant function L : TG x V* — > R relative to the 
affine action 

(v g , a) ^ (/iw 9 , fea + c (/i)). (2.23) 

Define the Lagrangians 

L ao : TG -> R, I : T(G xfxV)^!, and L y : T(G xTxF)^!, 

(see (2.3), (2.11) ; (2.15)) and the reduced Lagrangians 

I : g x V* -> R and / v : x T(V* xV)^l 

(see (2.4), (2.17)). Let (a(t),n(t)) and (a (£), ^o(^)) ^ e curves in V* x V related by the 
condition 

iP(g(t),a (t),v (t)) = (g(t),a(t),v(t)). 
For simplicity we assume that g(0) = e. Then the following assertions are equivalent. 

i One of the four conditions of Theorem 2.1 ( affine Lagrangian semidirect product theory) 
holds, with a(0) = ao- 

ii The curve (g(t), ao(£), Vo(t)) is a solution of the Euler- Lagrange equations associated to 
L, with a (0) = ao (this implies a (t) = a is a constant). 

iii The curve (g(t),ao(t),vo(t)) is a critical point of the action 

h _ 

L{g,g, a ,v a ,a ,va)dt 

1 

for variations 5g(t),5ao(t), and Svo(t) vanishing at the endpoints. 

iv The curve (g(t),a(t),v(t)) is a solution of the Euler- Lagrange equations associated to 
L v , with a(0) = a . 



t, 
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v The curve (g(t),a(t),v(t)) is a critical point of the action 

L v (g,g,a,v, a,v)dt 

to 

for variations 5g(t),Sa(t), and Sv(t) vanishing at the endpoints. 

vi The curve (£(t),a(t),v(t)) is solution of the Lagrange-Poincare (or reduced Euler- 
Lagrange) equations (2.10) associated to l v , with a(0) = a®. In our case, these equa- 
tions read 

d 5l v f5l v \ 

dtl* (2 24) 

dt 8a 5a ' dt 8v 5v 
vii The curve (£(t),a(t),v(t)) is a critical point of the action 

h 

l v (£, a, v, a, v)dt 

to 

for variations 5a(t), and 5v(t) vanishing at the endpoints, and variations <5£ of the form 

where n(t) e g vanishes at the endpoints. 

Proof. This theorem is a direct consequence of the discussion above and of the Lagrangian 
reduction theorem (Theorem 2.2). To see this, it suffices to show that the constrained 
variational principle in vii coincides with the one in Theorem 2.2 iii for the particular case 
when Q = G x V* x V. Indeed, we can identify the equivalence class [(g, a, v), £]g £ = 
(G xV* xV) x G g with the element (a, v, Ad g -i £) G V* x V x g. Using the trivial connection 
A(g, g, a, v, a, v) = gg~ x , the associated covariant derivative in q = V* x V x g is readily seen 
to be given by 

^- t (a(t),v(t),at)) = (a(t),v(t),at)). 

Thus, since the curvature of A vanishes, the general expression 5 A v = jyrfj + [v, fj] + B(Sx, x) 
reads simply 

This proves the result. ■ 

Remark 2.4 Note that the curves v(t), v Q (t) do not appear in point i. However, we can 
always recover them. The equation for v (t) is given by (2.14). The evolution for v(t) is 
implicitly contained in the definition of ip, namely, v(t) = g(t)~ l vo(t), where Vo(t) is the 
solution of (2.14) and g(t) is the G-part of the base integral curve for the Lagrangians L ao , 
L v , or L (they all give the same g(t) by construction). Equivalently, the equation for v(t) 
can be obtained by writing the Euler-Lagrange equation for L v relative to the variable a. 
An easy computation yields 

dL 

v - g' 1 -^ + (g~ l g)v = 0. (2.25) 
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We now check directly the equivalence between parts vi and vii. If 



S := l v (^,a,d,v,v)dt = / (l(£,a) + (v,a + £a + dc(£))) dt 

J tr, J to 



we get 



ti 
t 



*S= I [(^,Sn + ^,Sa\ + (5v,d + ^a + dc(0) 

+ (v, 5a) + (v, (5£)a + £(5a) + dc(5f)> ) dt 



ti 

to 



SI \ /si 

5£) + (—,5a) + (5v, d + £a + dc(£)> 



to 



<5£ / \5a 

- (v, 5a) - (£v, 5a) - (v o a, 5Q + (dc T (v), 5£)j dt 

~d\s\ + ad ^ Ji + dl^ v ° a ~ dcT ^^ - a ^l( voa - dc T 0)), v 

+ (5v, d + £a + dc(£)) + — v — £u, <it, 

where, in the last equality, we have used the relation 5^ = f] — [£,rj\. Thus, 5S = is 
equivalent to 

51 

d + £a + dc(f) = 0, — - v - £v = 0, 

da 

and ^ 

— — = ad, — + — (u o a - dc T (w)) - adUv o a - dc T (v)). 
dtd£ 5 d£ dt 5 

Using (2.22) and the the first two equations, the preceding equality reads 

d 51 ,„ 5/ 5/ , T (51 

-TT7 = ad f 77 + T - ° a - dc T - 

dr 04 ^4 oa \oa 

Thus, the variational principle in vii holds if and only if the following system 



d 51 ^51 51 , T ( 51 

— — = adt — + — o a - dc — 
dt5£ ? 5£ 5a \5a 

SI 

a + £a + dc(f) = 0, v - £v = 

oa 



(2.26) 



is satisfied. These equations are readily seen to be the Lagrange-Poincare equations associ- 
ated to I , and are equivalent to the affine Euler-Poincare equations (2.8), together with the 
advection equation for the variable a, and the equation for v. 

Remark 2.5 Recall the formula 
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If the horizontal equations for a and v hold (the second line in the system (2.26)), a direct 
computation shows that we have the equivalences 

5l v d 51 „5l 51 , T f5l 

— — = constant — — = ad f — + — oa-dc — 

d£ at o£ ? 04 5a \5a 



d 5l v 5l v 
51 



v 



5jv_ 
St 



where g(t) is determined by g{t) g(t) = £(£), g(0) = e. The last equivalence is a direct 
consequence of the general formula 

j t (Ad* Htyi p(t)) = Ad* Htyi (JLp(t) - sd* m n(t)j , m = Hty'hit), 

where t 1— > h(t) and t *—>■ p(t) are smooth curves in G and g*, respectively. 

Let p G g* and & £ 7* be given. Assume that the linear system w o b — dc T (w) = p 
has always a solution. Under these circumstances we can always choose the initial condition 
v(0) G V such that the constant in 5l v /5£ = constant vanishes. This is a convenient initial 
condition because it is preserved by the flow. One can think of this initial condition giving 
rise to 5l v /5C, = geometrically as one does with the vanishing angular momentum condition 
in the falling cat problem. 



3 Continuum spin system 

In this section we apply the previous theory to the dynamics of spin systems. 

The order parameter for a spin system is any Lie algebra, such as = su(2). Thus, 
besides the usual continuum degrees of freedom, Hamilton's principle for a spin system will 
depend on smooth functions of space and time that take values in a Lie algebra 0. Obvious 
examples are the magnetization vectors of magnetized fluids and spin-glass fluids, and a case 
can be made for also regarding He 3 -A and He 3 -B as spin systems. Another example, would 
be a fluid of oriented nano-particles, whose order parameter would take values on the unit 
sphere. Yet another example would comprise a fluid system of particles, each of which is a 
2- level qubit. The order-parameter dependence of the spin-system's prescribed interaction 
energy determines its Hamilton's principle. 

It is shown in Gay-Balmaz and Ratiu [2009] that the equations of motion, as well as 
their variational and Hamiltonian structures, can be obtained by affine Euler-Poincare or 
Lie-Poisson reduction. Given a manifold T> and a Lie group O, we consider the Frechet Lie 
group T(V,0) of smooth functions x '■ D ~^ ®i an d the vector space V* = fi 1 (P, 0) of 
one-forms with values in the Lie algebra of O. When T> is endowed with a volume form p,, 
this space is in a natural way the dual of the space V = X(V,0*) of 0*- valued vector fields 
on D, the duality pairing being given by contraction and integration over T>. Let \ G G act 
on 7 G V* via the left affine representation 

7 h-> Ad x 7 + x^X" 1 - 
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This action is natural in the sense that it coincides with the usual gauge transformation of 
a connection on the trivial principal bundle O x T>. The Lagrangian of a spin system is a 
function 

L^ \ TG — ► R, 

such that L(z/0,7) := L y (u^) is invariant under the left affine action of \ G G on TG x V* 
given by 

0i/» 7) >-> (x*ty, Ad x 7 + xdx _1 ) • 

We can now apply the general theory of affine Euler-Poincare reduction to the group G = 
J-'iV, O), the space V* = 0), and the group one-cocycle c G ^(G, V*) given by 

c(x) = X^X^ 1 - 
One can check that the cocycle property 

c(xV0 = c(x) + Ad x c(^) 

holds, relative to the representation 7 1— > Ad,/, 7 of G on 7*. As in the abstract theory, we 
use the notation # x 7 = Ad x 7 + x^X 1 f° r the affine representation. 
For w G V and 76 V*, the associated diamond operation is given by 

w o 7 = Tr(ad* u>) = ad*, w' 1 . 

Since c(%) = x&X 1 for ^ G and w £ V, from (2.2) we obtain the relations 

dc(u) = —dv and dc T (w) = divw. 

Thus, the expression for the cotangent-lift momentum map (2.19), 

d c T ( — \ - — o - J ( — 

appearing with the opposite sign in the right hand side of the affine Euler-Poincare equations 
(2.8), in this case reads 

^-*(<)=4S- (3l) 

Remarkably, the cotangent- lift momentum map in this case recovers an expression known as 
the covariant divergence associated to 7, defined otiwG X(TJ, 0*) by 

div 7 w = div w — Tr (ad* w) . (3.2) 

Making use of these observations, the affine Euler-Poincare equations (2.8) become 

— — - ad — + div 7 — = 0. 3.3 
at bu ov 07 

Since the infinitesimal action of v on 7 induced by the representation 7 i— > Ad x 7 is 

~{v = ad v 7, 
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the advection equation (2.5) for 7 is 

J^ T -dx/ + ad„7 = 0. (3.4) 

Here also, the expressions associated to the affine representation combine in a nice way, since 
the quantity (du — ad„ 7) is known as the covariant derivative of v associated to 7, and is 
usually denoted by d 7 zA The 7-advection equation (3.4) can thus be rewritten as 

| 7 = cTV. (3.5) 

In the particular case V = M 3 and O = 5*0(3), equations (3.3) and (3.5) appear in the context 
of the macroscopic description of spin glasses, see equations (28) and (29) in Dzyaloshinskii 
[1980] and references therein. See also equations (3.9), (3.10) in Isaev, Kovalevskii, and Peletminskii 
[1994], system (1) in Ivanchenko [2000] and references therein for an application to magnetic 
media. In this context, the quantity 51 /5v is interpreted as the spin density, v is the in- 
finitesimal spin rotation and the curvature B := d 7 7 is the disclination density. 

Remark 3.1 Thus, the covariant derivatives appearing in the Hamiltonian operator as gen- 
eralized 2-cocycles in the dynamics of spin systems arise directly from the affine Lie symmetry 
of the Hamilton's principle. 



Variational principles for spin systems. We now specialize to the case of spin systems 
some results obtained abstractly in Theorem 2.3. Start with a Lagrangian 

L 70 : TT{V, O) -> K, 

describing the dynamics of a spin system and depending on a parameter 70 G V* = Vt l (T>, 0). 
Suppose that the function L, defined by L[y^, 70) := L l0 {v^), is invariant under the left affine 
action of \ G O) on TTp, O) x tt\V, 0) given by 

(z/0, 7) ^ (xi/0, Ad x 7 + xd^T 1 ) • 

We will also use the notation instead of u^. The Lagrangians L and L v (see (2.11), 

(2.15)) defined on T {Tip, O) x Vt l (V, 0) x X(V, 0*)) are given by 



L(^,^, To, 70,^0,^0) = L(V>,^,7o) + / (7o -wo) fi 

and 



v 



L v {i),i),^, A i,w,w) = L(ip,ip,9^) + J wi^i-d 1 ^ V)/^- 
The reduced Lagrangian l v (see (2.17)) defined on Tip, 0) x T(fi 1 (P, 0) x X(T>, 0*)) is 

l v (y, 7, w,j, w) = l{v, 7) + / w-(j - d 7 z/) fi. 



v 



We are now ready to state several variational formulations of the dynamics of spin systems 
obtained from Theorem 2.2. 



Gay-Balmaz, Holm Sz Ratiu Variational principles, spin systems and Kir chhoff rods 16 



Let x(t) be a curve in the group G = J-(T>, O) and fix an element 70 G V* = Q}{V, 0). The 
choice 70 = is allowed and important. Define the curve u{t) := G T(T>, 0). Let 

(7o(^), u>o(t)) and (j(t),w(t)) be two curves in Q l (V, 0) x X(V, 0*) related by the conditions 

Md*)- 1 7o(*) + X(*) _1 dx(*) = 7(*) and Ad x(t) ^o(*) = w(t), 

and such that 70 (0) = 70. Suppose, for simplicity, that x(0) = e. Then the following are 
equivalent: 

i Hamilton's variational principle 

s / L 10 (x,x)dt = 0, 

holds, for variations Sx{t) vanishing at the endpoints. 
ii Hamilton's variational principle 

<5 J \L(x,X,lo) + J io-w^jdt = Q, 

holds, for variations Sx(t), Sjo(t), and 5wo(t) vanishing at the endpoints. 
iii Hamilton's variational principle 



5 



J t 1 (l (x, X, d x l) + jf w- (7 - cr*" 1 *)) dt = 0, 



holds, for variations 5x(t),^7(t), and vanishing at the endpoints. 

iv The reduced variational principle 

5 [ 1 (l(u, 7) + / (7 - <Tv)\ dt = 0, 



Jt \ Jv J 

holds, for variations Sj(t), and 5w(t) vanishing at the endpoints, and variations 5u{t) 
of the form 

dr/ 

where rj(t) vanishes at the endpoints. 
v The reduced variational principle 

5 [ l(v, 7)df = 0, 
holds, for variations 8u(t) and 5~f(t) of the form 



dri 

8v = — — h [u, 77] and £7 = dCn 



where n(t) vanishes at the endpoints. 
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Note that the curve 7(f) is not present in part i, but it can be recovered from the curve 
x(t), through the relation 

7(f) = x(t) -i7o = Ad x(t) -i 70 + x{tY l d X {t). (3.6) 

The curve 70 (t) = 70 is a constant. The curve w(t) is not present in the parts i and v, but 
it can be recovered by solving the equation 

w — ad w — — = 0, 

07 

obtained from stationarity under variations with respect to 7. The variable w is the momen- 
tum canonically conjugate to the 1-form 7, just as the electric field is canonically conjugate 
to the vector potential in electromagnetism and its Yang-Mills generalization for any gauge 
group. The w-equation may be written equivalently as 

w - ad H w - Ad —— = . 

O'-f 

The curve Wo(t) is related to w(t) through the relation 

w (f) = Ad* (t) _i w(t), 
or it can be obtained directly by solving the equation 

w - — = 0. 

o~fo 

All these variational principles are equivalent to the afflne Euler-Poincare equation (3.3) 
together with the equations for 7 and w. 
The conservation law (2.9) for spin reads 

9 < Ad- , «) + div- ( x «) = 0, (3.7) 



dt \ x 8v J V ^7 / 

where 70 is the initial value of 7. In the important case when 70 = 0, one finds the conser- 
vation law, 

IKS 

cf. §10 in Holm [2008], especially equation (10.1.6). 

Remark 3.2 The curvature of 7(f) corresponds to the density of defects in the order pa- 
rameter. If initially 70 = 0, then the curvature of 7(f) remains zero for all t (see (3.6)) and 
no defects arise. However, a non-zero initial value for the connection 70 7^ allows nonzero 
curvature, corresponding to defects, which in turn alter the spin dynamics of the system 
(3.7), now rewritten equivalently as, 

^d-^+div^Wfad. ( x «)). (3.8) 



dt V x <W V h J V 70 V h 

The right hand side of this equation shows the effect of a non-zero initial value of the 
connection 7 on the dynamics of the spin density. 
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Remark 3.3 Using the advection equation (3.4), that is, ^ = du + [7, u], if c is a closed 
loop in V, we get 

~d~i j[ 7 = ^ dl/ ~ adu 7) = - j^ad„ 7. 
In particular, if O is Abelian, we get 

Returning to a general Lie group O, we calculate the dynamical equation for the defect 
density B := d 7 7 = d7+ [7,7] by taking the 7-covariant exterior derivative (covariant curl). 
We get 

B = d7 + [7, 7] + [7, 7] = d (du + [7, u\) + [dp + [7, u], 7] + [7, du + [7, 1/]] 
= [d 7 + [ 7 , 7 ],z/] = [£,!/]. 

Thus the advection equation for the curvature reads 

dB 

— + ad u B = 
at 



which implies 

d 
dt 



s 



B = - &d u B, (3.9) 



for any compact two-dimensional submanifold S of T>. 

In particular, if O is Abelian, the advection equation for the curvature becomes the 
conservation law dB/dt = 0. For example, if dimP = 3 and O = S 1 , we can define the 
vector field B = (*!?)" and the integral of the left hand side of the previous equation is the 
flux of B; we get in this case 

|// s (b.„ms = o. 

For a general Lie group, equations (3.8) and (3.9) govern the nonlinear coupling between 
spin and defect density dynamics. If initially 70 = then (3.6) shows that the curvature 
B(t) of 7(t) remains zero for all t and no defects arise. However, a non-zero initial value 
for the connection 70 7^ allows nonzero curvature, corresponding to the initial presence of 
defects. The defects alter the spin dynamics of the system in (3.8), which in turn feed back 
nonlinearly to influence the defect dynamics in (3.9). The right hand side of equation (3.9) 
shows the non-commutative effect of the spin density on the defect density due to a non-zero 
initial value 70 7^ of the connection 7. Namely, non-commutativity J J ad„ B^O generates 
defects when 70 7^ 0. 



4 Elastic filament dynamics and Kirchhoff 's theory 

The equations of motion for continuum mechanical systems can be formulated in three 
representations: material or Lagrangian, spatial or Eulerian, and convective or body. For 
filaments, there is a fourth representation, due to Kirchhoff, that is often very convenient. 
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4.1 Kirchhoff 's theory and spatial representation 

In this subsection we briefly recall the relevant aspects of Kirchhoff's theory of rods (see 

Dichmann, Li, and Maddocks [1996] for a modern formulation). We shall apply it here to 

models that depend explicitly on order-parameter variables. In Ellis, Gay-Balmaz, Holm, Putkaradze, and 

[2009] one can also deal with models that have nonlocal interactions by using affine Euler- 

Poincare reduction instead of Kirchhoff's theory. 

The Kirchhoff rod is described by a vector r(s) and an orthonormal basis of director 
{di(s), d 2 (s), d 3 (s)}. The curve r(s) is interpreted as the configuration of the center line, 
the triad {d^} can be interpreted as providing information concerning the orientation of the 
material cross-section of the rod. Given a fixed orthonormal basis {Ei,E2,E3}, we define 
the transformation matrix A(s) G 5*0(3) by requiring the equality 

d<(s) = A*(s)E fcj z = l,2,3. (4.1) 

From the variables r and A we define the body angular velocity uj and the linear velocity 7 
by 

Q = A -1 A and 7 = A" 1 !-, (4.2) 

where : (IR 3 , x) — > (so (3), [,]) is the Lie algebra isomorphism given by uv = u x v for all 
v G M 3 . Thus, in an orthonormal basis of M 3 and u G M 3 , the 3x3 antisymmetric matrix 
u := u G so (3) has entries 

u jk = (u) jk = -e jk iu l . (4.3) 
Note that the relations (4.2) can be rewritten using the group structure of SE(3) as 

(^, 7 ) = (A,r)" 1 (A,f). 

The linear momentum density p is defined as p(s) = pd(s)'y(s), where pd(s) is the local mass 
density of the rod. In that case, the kinetic energy due to linear motion Ki in is given by 

Kun = \j Pd (s)\\r(s)\\ 2 ds = ~J AMUA-^fda = ~ J p d (s)\\~r(s)\\ 2 ds. 

Consequently, we have 



SKu 



I'm 



5j 

The local angular momentum in the body frame (di(s), d2(s), da(s)} is defined by 7r(s) : = 
I(s)w(s), where I(s) is the local value of the inertia tensor expressed in body coordinates. 
Relative to body coordinates we have 7r*(s) := ¥j(s)u : '(s). Thus the local kinetic energy due 
to rotation is given by 

K rot = \ J w(s)-I(s)wds 

and hence 

tt 5K rot 

7T = = — . 

0<*> 

In Kirchhoff's approach, the conservation laws are written in terms of variables expressed in 
the fixed spatial frame {Ei,E 2 ,E 3 }. To distinguish it from p and iz which were expressed 
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in the body frame (di(s), d2(s), da(s)}, we shall denote the same vectors in the fixed spatial 
frame {Ei,E 2 ,E 3 } by p (E) and 7r (E) . Thus, (4.1) yields the relations 



7T = 7r l di = 7r (E) ' fc E fc and p = p% 



P 



(E),fc E 



where 



A 



SK, 



rot 



and p (E) ' fc = p d A*Y 



A 



5K, 



i in 



5j 



(4.4) 



Thus, the vectors p^ E ^ and of body linear and angular momenta expressed in the spatial 
frame are connected to the local body quantities as 



tt^ = A 



5K 



rot 



6u) 



and p (E) = A 



6K, 



l in 



5~f 



(4.5) 



In general, it is assumed for physical reasons, that the Lagrangian in Kirchhoff's formu- 
lation has the form 



7, ^, r) = K lin (7) + K rot (u) - E(n, r) , 



(4.6) 



where the potential energy E(fl, T) is usually taken to be a quadratic function of the defor- 
mations Q := A _1 (9 S A and T := A~ 1 d s r, but more complex expressions are possible as well; 
we shall not restrict the functional form of that dependence. In this case, the body forces 
n = SI/5T and torques m = 51/ 5Q are connected to the transformed quantities n^ E ^ , 
in Kirchhoff's theory as 

SI ™ . 51 



li 



( E ) = A—, m( E ) = A—. 



(4.7) 



The balances of linear and angular momenta in Kirchhoff's approach (cf. eqs. (2.3.5) and 
(2.3.6) of Dichmann, Li, and Maddocks [1996]) are given by 



<9p( E > dnW _ 

dt ds 
dirW <W E > dr 



ds 



Ot 



ds 



x n 



(E) 



1. 



(4.8) 
(4.9) 



where f and 1 denote the external forces and torques, respectively. 

We now rewrite equations (4.8), (4.9) in terms of spatial variables, using the formula 



Ad( A , r) -! (/x, rj) = (A/x + r x Ary, Arj) . 



(4.10) 



for the coadjoint action of SE(3) on se(3)*. The spatial momenta - denoted by a superscript 
(S) - become 



(TT^p^Ad^j 



/ 51 51 
\ 5oj ' 5~f 
(Tr^ + rxp^.p^l) 



51 51 . 51 

A— + r x A— , A— 
5u 07 07 



(4.1i; 
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and the spatial torques and forces are 

(n «.,„«.) = = AlW , 1 («,«) = (A« + rxA«, A «) 
= (m( E )+rxn( E ),n( E )). 
The conservation laws in Kirchhoff theory may now be written as 

|(irW,pW) + !(mW,n«>) = (T,r), 

where T := r x f + 1. 

4.2 AfRne Euler-Poincare formulation 

We now show that, in the case of potential forces, we can obtain these equations by affine 
Euler-Poincare reduction. 

Consider the group G = J-"([0, L], SE(3)) of all smooth maps from the interval [0, L] 
into the group SE(3), whose elements are denoted by (A, r). The space of affine advected 
quantities is chosen to be the dual vector space 

V* := fi 1 ([0,L],se(3)) ®^([0,L],M 3 ) 3 (fi,I\p). 

Note that since [0, L] is one dimensional, the space f2 1 ([0, L],se(3)) can be identified with the 
space J-"([0, L],se(3)) of smooth maps with values in se(3). We choose to think of these func- 
tions as one-forms, since it is that interpretation that generalizes to the higher dimensional 
case. We consider the affine representation 9 of the group G on V*, given by 

( a,,)(Q, r, p) = (A, r )(n, r, p) + C (A, r), 

where the first term denotes the representation of G on V* defined by 

(A,r)(n,r,p) = (Ad (A , r) (n,r),Ap) (4.14) 

and the second term is the group one-cocycle given by 

c(A,r) = ((A,r)9.(A,r)- 1 ,-r), (4.15) 

where s 6 [0,L]. The fact that c is a group one-cocycle relative to the representation 
(4.14) (that is, it satisfies (2.1)) follows from a direct computation. Note that the first two 
components of c form the left version of the cocycle appearing in the theory of complex 
fluids, see Gay-Balmaz and Ratiu [2009]. The infinitesimal linear action of the Lie algebra 
element (u>,7) G JF([0, L], se(3)) is 

(w, 7)(fi, T, p) = (ad( W)T )(n, r),wp) = (w x f2, a; x T - fi x 7, w x p) 

and the associated diamond operation is 

(u, iu, /) (fi, r, p) = (u x fj + w x r + / x p, iu x n). 



(4.12) 
(4.13) 
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Concerning the cocycle c, we have the formulas 

dc(w, 7) = (-d s u, -dgj, -7), dc T (u, w, f) = (d s u, d s w - f). 
Using all these expressions, the affine Euler-Poincare equations (2.8) become in this case 



, 51 , , , 51 51 51 51 

(d t + wx — + (d s + fix — = — x 7 + — x T + — x p 
5u> dil 07 ol dp 

\ <^ /n \ ^ <^ 

O i + W x)- + ( a i + nx)- = -. 

The advection equation (2.5) for Kirchhoff's theory are 

ti + u> x ri = 

f + w x r = 9 S 7 + o x 7 



(4.16) 



(4.17) 



If we fix the initial values Qq, Tq, p , the conservation law (2.9) reads 



0_ 

dt 



51 51 
(A ' r)_1 V^' ^7 



Ad*A ^-1 



div (Oo,r ) 



Ad*A 



51 51 



(A ' r)_1 Un'ar 



or, more explicitly, 



Ad* ^ 



ds 



51 51 



A M + rxA ^J X ^ + A ^ Xro + A ^ X ^ A ^ Xr ° + A ^J 
= (m^ x ft + x r + A^ x p , x T + A^ =: (T, f ) (4.18) 

using (4.12). Note that the right hand side of this equation is formed by an external torque 
T and an external force f that recover (4.13) in terms of the spatial variables 

=^ («.,«) and (m«,n«)=AdJ M .,g,^). 



If f2 = 0, r = 0, p = 0, the conservation law (4.18) simplifies to 



d_ 
di 



Ad; 



51 51 



d_ 

ds 



Ad; 



51 51 



0,A 



5l_ 

Tp 



(4.19) 



Note that in this case, there is no external torque. Conversely, non-zero initial values Qq, 
r , p Q produce the external torque in (4.18). This torque phenomenon for the Kirchhoff 
rod corresponds to the effect of non-zero curvature (defect density) in Remark 3.2 for spin 
systems. The corresponding connection for the Kirchhoff rod is given by the pair (fl, T). 
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When the derivative 8l/8~f is proportional to 7, which is the case for Kirchhoff 's La- 
grangian, equations (4.19) can be rewritten as 

at as dp 

9 (E) 9m (E) 3 « 

— + — + 9j Xn,, = A ^ X1 -' < 4 - 21 > 

and we recover the balances of linear and angular momenta in Kirchhoff's approach when 
central potential forces are considered; that is, when the potential depends only on the 
magnitude of distance \r\ = \p\. The variable p = A~ 1 r is the vector r from the center of 
coordinates to a point on the rod, as seen from an orthogonal frame fixed at that point on 
the moving rod. 

Returning to the general case, recall that the reduced Lagrangian is of the form 

z(«, 7, ft, r, P ) = K rot (u) + # Kn ( 7 ) + E{n, r, P ). (4.22) 

Since \r\ = \p\, one can replace a central force in r by one in p. For example, one could 
discuss the dynamics of a Kirchhoff rod moving in a gravitational field, such as a tether 
hanging from an orbiting satellite and feeling the force of Earth's gravity. 



Generalizations. In fact, our approach generalizes, without additional difficulties, to the 
case when the variables A and r are defined on an arbitrary manifold T> and take values in 
an arbitrary Lie group O and an arbitrary vector space E, respectively. Correspondingly, 
the other variables are 

fieO^o), reQ^E), P eF(v,E), 

where o denotes the Lie algebra of O. The symmetry group G = ^([0, L], SE(3)) of Kirch- 
hoff's theory generalizes to the group G = J 7 (T>,S), where S is the semidirect product 
S = 0®E relative to a left representation of O on E. Thus, the group multiplication is 
given by 

(X, v) (ip, w) = (xV>> v + xw) ■ 

The space of affine advected quantities is chosen to be the vector space V* := Q l (V,s) © 
F(V,E) 3 (fl,r,p), dual to the space V = X(V,s*) © T{V,E*) 3 (m,nj). The affine 
representation 9 of the group G on V* has the same expression as before, namely, it is given 
by 

(A , r) (fi,r» = (A,r)(fi,r,p) + c(A,r), 
where the first term denotes the representation of G on V* defined by 

(A, r )(n, r, p) = (Ad (A|r) (n, r), a p ) = (Ad A n, at - (Ad A n)r, a p ) 

and the second term is the generalization of the cocycle (4.15) given by 

c(A, r) = ((A, r)d(A, r)~\ -r) = (AdA" 1 , -dr - (AdA^r, -r). 
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Using the formulas 

(u, w, f) o (fi, T, p) = (ad^. u* + w i o Ti + f o p, -Q^) 

and 

dc(u, 7) = (-da;, -d7, -7), dc T (w, to, /) = (div(tt), div(w) - /), 
the affine Euler-Poincare equations read 

9 \ 5/ o SI 51 
m +iU )^ + dw ST = v 



(4.23) 



where div is the covariant divergence defined in (3.2). The advection equations (2.5) become 
in our case 

n = d n uj , 

f + uT = d n 7 , (4.24) 
^ p + ^p = 7, 

where 

d u; = did + ad^ a; and d^7 = d7 + Q^y. 

When specialized to the case V = [0, L] and S = SE(3), these equations recover (4.16) and 
(4.17) of Kirchhoff 's theory. 

4.3 Variational principles for the Kirchhoff rod 

In this paragraph, we specialize some results obtained abstractly in Theorem 2.3 to the case 
of the Kirchhoff rod. In order to give a more transparent vision of the underlying geometric 
structures, we consider the generalization described above, that is, we replace the interval 
[0, L] by an arbitrary manifold T> and we replace SE(3) by the semidirect product S = O (S) E 
of a Lie group O with a left representation space E. 
Start with a Lagrangian 

L(n ,r ,p ) '■ TTiV, S) -> R, 

depending on the parameters 

(fi ,r ,p ) e n\v, s )®T(v,E). 

Suppose that the function L, defined by 

L(A, A, r, r, fi„, T , Po) := £(n ,r 0l/ 9o)( A > A > r > 0» 
is invariant under the left affine action of (%, u) G 5) on 

(A, A, r, r, fi, I» G TF{V,S) x (^(D, s) © ^(P, £)) 
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given by 



" A " 




yA 


A 
i\ 




vA 


r 




v + x r 


r 




Xr 


n 




Ad x n + yd* -1 


r 




X V-(Ad x Q)v 


. p . 




_XP-v 



This seemingly complicated action is nothing else than the affine action (2.23) written for our 
G := T(V, S) and V* := Q 1 (V , s) © T '{V , E). Such a function L is completely determined by 
the expression of the reduced Lagrangian /, for example the Lagrangian (4.22) of Kirchhoff 's 
rod. 

The Lagrangians L and L v defined on the tangent bundle 

T[F(V,S)x (Q\V,s)®^{V,E)) x (X{V,s*) ®F(V,E*))] 

are given by (see (2.11), (2.15)) 

L(A, A, r, r, Q , T , p , tt , t , p , u , w , f , u , w , f ) 

= L(A, A, r, r, Q , r , po) + / (l -u + t -w + po-fo 

Jv 

and 

z7(A,A, r, f, q, r, p, n, r, p, u, w, f, if, w, f) 

=L{A, A, r, f, Ad A Q + AdA^ 1 , Af - (Ad A Q)r - d AdA ~V, Ap - r) 



v 



u-[n-d u u) + w-{T + LuT-d u -f) +f.(j, + U p->y) 



where u = A X A and 7 = A l f. The reduced Lagrangian l v defined on 

F{V,s) x T [(^(P,*) @F{V,E)) x (£(£>, s*) ®.F(D,£*))] 

is (see (2.17)) 

Z v (a;,7, fi, T, p, fi, f, p, it, iy, /, it, tu, /) 



=i(w,7,n,r,p) 



V 



u-(£l-d n u) +w(t + LuT-d n -f) + f-(p + up-j) . 



We are now ready to state the various variational formulations of the dynamics of Kirchhoff 's 
rod obtained from Theorem 2.3. 

Let (A(t),r(t)) be a curve in the group J-'iV^S) and fix an element (Q ,T ,p ) £ 
r2 1 (D,s) © TiT),E\ The choice (f2 ,r ,po) = (0,0,0) is allowed and important. Define 
the curve 

(w(*),7(<)) := (A(t)^A(t), A(t)-V(t)) £ ;F(Z>,s). 
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Let (Q (t),T (t),p (t),uo(t),w (t), f (t)) and (Q(t),T(t), p(t),u(t),w(t), f(t)) be two curves 
related by the conditions 

(n, T, p) = (A, r) _1 (fi , r , p ) and {% w, f) = (A, r)'\u , w , f ) 

and such that 

n (o) = n , r (o) = r , p (o) = p - 

Suppose, for simplicity, that (A(0),r(0)) = (e,0). Then, the following are equivalent, 
i Hamilton's variational principle 

S I ^(n ,r ,p ) ( A > A > r ' r)dt = 0, 
J t 

holds, for variations 5A(t),Sr(t) vanishing at the endpoints. 
ii Hamilton's variational principle 



to 



L{A, A, r, f, fi , r , p ) + / (lo-u + t -w + po-fo) dt = 



D 



holds, for variations 5A(t), Sro(t), 5Qo(t), 5To{t), 5po(t), 5uo(t), 5wo(t), and 5fo(t) van- 
ishing at the endpoints. 

iii Hamilton's variational principle 

5 J 1 [l (A, A, r, r, Ad A Q + AdA" 1 , Af - (Ad A fi)r - d AdA ~V, Ap - r 



v 



U-[Q-d n Lj) +W[t + LUT - d n 7 ) + f.(p + U p-y) 



where uj = A 1 A and 7 = A 1 r, holds, for variations SA(t), Sr(t), 5Q(t), ST(t), Sp(t), 
Su(t), 5w(t), and Sf(t) vanishing at the endpoints. 

iv The reduced variational principle 



to 



s / i(uj,y,n,r,p) 



+ 



V 



u-(n-d n u) + w(t + uT~d n y) + f-(p + up-y) 



dt = 0, 



holds, for variations 5Q(t),5T(t),5p(t),5u(t),Sw(t), and 5f(t) vanishing at the end- 
points, and variations 8u(t) and Sy(t) of the form 

drj dv 
5uo = — + [u, n] and Sy — — + uv — 777 



where rj(t),v(t) vanish at the endpoints. 
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v The reduced variational principle 



5 / 1(uj, 7, Q, r, p)dt — 0, 
Jt 



holds, for variations 5u{t),S'j(t),SQ(t),ST(t),5p(t) of the form 

r ^rj dv 

OU) = — + [U, T}\, 5 7 = — + UV - 777, 

8fl = d Q rj, 5T = d Q v — r/T and 5p = v — np, 
where n(t),v(t) vanish at the endpoints. 

Note that the curves fl(t),T(t), p(t) are not present in part i, but they can be recovered 
from the curves A(t),r(t), through the relation 

(£l(t),T(t),p(t)) = 0(A(t) )r (t))-i(no,r o ,/0o). 

The curves Qo(t) = Qo,T (t) = r ,po(^) = Po are constant. The curves u(t),w(t),f(t) are 
not present in parts i and v, but they can be recovered by solving the equations 

u-ad*u + 7 ou;--^ = 0, 

61 n 
w + uw - — = 0, 

ol 

51 



f + uf-- = 0, 
8p 



or, equivalently, the equations 



u - adl_ 1A u + (A-V) o w - Ad^ (H - r o = 0, 

BL 

w + (A- 1 A)w-A' 1 — = 0, 

OL o 

/ + (A- 1 A)/-A- 1 — = 0. 

dp 

The curve (uo(t), Wo(t), fo(t)) is related to (u(t),w(t), fit)) through the relation 

(u ,w ,f ) = (A, r)(u,w,f) 
or it can be obtained directly by solving the equations 

U ~ "nrT = °' W 0~ = and JO - 7^— = 0. 



All of these variational principles are equivalent to the affine Euler-Poincare equation (4.23) 
together with the equations for Q, T, and p. 



Gay-Balmaz, Holm Sz Ratiu Variational principles, spin systems and Kirchhoff rods 28 



The case of Kirchhoff' s rod. We now rewrite the points iv and v in the particular 
case T> = [0, L\ and S = SE(3), that is, the case of Kirchhoff theory. We thus find that 
Kirchhoff 's equations for a Lagrangian I are equivalent to the following. 

iv The reduced variational principle 

u- ( S7 — 8 h uj — Q x uj 



5 jf 1 ^(w, 7 ,n,r,p) + jJu-[ 

+w- (t + uj x r - d s j - ft x 7^ + /(p + uj x p - 7) J eft = 



holds, for variations 8fl, 5T, 5p, 8u, 8w, and 5f vanishing at the endpoints, and 
variations 8uj and £7 of the form 

ouj = — — h uj x 77 and 07 = — — \- uj x v — 77x7 

where rf(t),v(t) vanish at the endpoints. 
v The reduced variational principle 



5 f 1 i(u,<y,n,r, P )dt = o, 

J t 



holds, for variations 5uj, 5~y, 5 ft, 8F, 5p of the form 

. dr) dv 

ouj = — — h uj x 77, 07 = — — \- uj x v — 77x7, 

= <9 S ?7 + 17 x 77, 5r = 9 s v + J7x7j + rx77 and bp = v — 77 x p, 
where rf(t),v(t) vanish at the endpoints. 

The variables (u, w, /) are conjugate momenta to the constraint variables (ft, T, p), respec- 
tively, in the reduced variational principle in point iv. 

Recall from the general theory developed in Section 2, that the integrand appearing in the 
action principle in iv above is the reduced Lagrangian I associated to the cotangent bundle 
Lagrangian L v , by standard Lagrangian reduction TQ — > TQ/G, see (2.15) and (2.17), for 
the definition of these Lagrangians. In particular, (2.15) shows how one can construct L v 
from the Lagrangian L = L(A, A, r, r) of the Kirchhoff rod. 
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